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Abstract 

We discuss an extension to the Chandrasekhar-Fermi method for the eval- 
uation of the mean magnetic field strength in molecular clouds to cases where 
the spatial orientation of the field is known. We apply the results to M17, using 
previously published data. 

Subject headings: ISM: clouds — ISM: individual (M17) — ISM: magnetic fields 
— polarization 

1. Introduction 

There exist few techniques that allow for the measurement of quantities that characterize 
the magnetic field in molecular clouds. At millimeter and submillimeter wavelengths, the 
orientation of the magnetic field is most commonly traced using polarimetry measurements 
from dust continuum emission (Hildebrand 1988). The strength of the magnetic field (in 
general, its line-of-sight component) can only be directly measured via the Zeeman effect 
(e.g., Crutcher et al. (1999); Brogan & Troland (2001)), usually at longer wavelengths. In 
order to gather as much information as possible about the magnetic field, the so-called 
Chandrasekhar Fermi (CF) method (Chandrasekhar & Fermi 1953) is often used to infer the 
strength of the plane-of-the-sky component of the field. Because this is achieved with the 
same polarimetry data that give the orientation of the sky-projected magnetic field, the CF 
method can act as a bridge between the polarimetry and Zeeman observations to provide an 
estimate for the magnitude of the mean field strength in a given cloud. 

In this paper, we discuss how a simple extension of the CF method can be used alone, 
i.e., without the need of Zeeman data, to infer the magnitude of the magnetic field; not 
only the strength of its plane-of-the-sky component. Furthermore, it will also be shown 
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that, contrary to the original CF method which only really works well when the magnetic 
field is located close enough to the plane of the sky, our generalization is valid regardless 
of the field's orientation in space. However, this can only be accomplished if and when the 
spatial orientation of the magnetic field is known. That is to say, not only the orientation of 
its projection on the plane of the sky is needed (from polarimetry), but also its inclination 
to the line of sight. This last piece of information can be obtained through the technique 
of Houde et al. (2002) which relies on the availability of spectroscopic measurements from 
suitable neutral and ionic molecular species, as well as polarimetry. 

Finally, we apply our extension to the CF method to already published data for the 
M17 molecular cloud (Houde et al. 2002), and infer a value for the magnitude of the mean 
magnetic field for this object. 



2. The CF equation 

It was originally asserted by Chandrasekhar & Fermi (1953) that the amount of disper- 
sion of the polarization angles measured from starlight (or dust continuum radiation) can 
reveal information about the magnitude of the magnetic field. With the assumption that 
the magnetic field is frozen to the ambient fluid, any (turbulent) motion within the gas in a 
direction perpendicular to the orientation of the magnetic field will be transmitted to, and 
distort, the field lines. Chandrasekhar & Fermi (1953) further assumed that such distur- 
bances would propagate as waves along the magnetic field lines at the Alfven speed, which 
they used as the starting point for their analysis. It follows that since dust grains are thought 
to be tied to the magnetic field lines (Mouschovias & Ciolek 1999), the amount of distortion 
in the field lines can be inferred from polarimetry. Similarly, the turbulent motion of the gas 
can be measured through the spectral line profiles of molecular species, for example. These 
two observed quantities are needed to evaluate the strength of the magnetic field through 
the CF method. 

Following, therefore, the original derivation of Chandrasekhar & Fermi (1953), we can 
write an equation for the mean value of the magnetic field B as 



B=V^, (1) 

where p and cr(v±) are, respectively, the mass density and the two-dimensional velocity 
dispersion (perpendicular to the field lines) of the matter coupled to the magnetic field. 
a (<p) is the dispersion in angular deviations of the field lines. Equation (1) is valid in the 
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small angle limit. 

In their estimation of the magnetic field strength in the spiral arms Chandrasekhar & Fermi 
(1953), identified a (0) with the dispersion in the orientation of the polarization vectors 
measured for distant background stars. Using the coordinate system of Figure 1 to define 
the spatial orientation of the magnetic field, with a the inclination angle of the field to the 
line of sight, and f3 the angle made by its projection on the plane of the sky, we find, for the 
case originally considered by Chandrasekhar & Fermi (1953), that 



a (0) = a (P) . (2) 

However, observations of this type probe only one direction in the lateral displacement of 
the magnetic field line. We must, therefore, make the following substitution for the velocity 
dispersion 



where a [v) is the total three-dimensional velocity dispersion of the gas (for cases of isotropic 
turbulence). Inserting equations (2) and (3) in equation (1) we obtain the original equation 
derived by Chandrasekhar & Fermi (1953) 



fl a(v) 

where B pos is the plane-of-the-sky component of the magnetic field (more on this below). 

Equation (4) is often used to measure the mean strength of the plane-of-the-sky com- 
ponent of the magnetic field in molecular clouds (e.g., Lai, Girart, & Crutcher (2003)). It 
has also been tested with magnetohydrodynamic (MHD) simulations to verify its domain 
of applicability (Ostriker, Stone, & Gammie 2001; Padoan et al. 2001; Heitsch et al. 2001; 
Kudoh & Basu 2003). Although the CF method has been found to work well for strong 
enough magnetic field, it also suffers from some shortcomings. Among these, is the fact that 
the equation (4) only really applies well when the magnetic field is located close enough to 
the plane of the sky. In fact, the method will fail when the field is aligned parallel to the 
line of sight (a = in Figure 1). 
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2.1. An extension to the CF method 

It would be desirable to extend the CF method to cases where the magnetic field is 
arbitrarily oriented in space. This, however, requires that observations can be made to 
measure not only (3 (the angle made by the projection of the magnetic field on the plane of 
the sky), but also a (the inclination angle of the field to the line of sight). Some methods 
have already been proposed to do such measurements. Myers & Goodman (1991) (see also 
Bourke & Goodman (2004)) modeled the magnetic field in molecular clouds with uniform 
and nonuniform components, and through a statistical analysis were able to evaluate the 
spatial orientation (i.e., they inferred a and f3) for the mean three-dimensional uniform field. 
More recently, Houde et al. (2002) have proposed a technique that combines polarimetry 
and ion-to-neutral line width ratio measurements (Houde et al. 2000a, b) to map the spatial 
orientation of the magnetic field across molecular clouds. This method has been used so far 
for three different objects: M17 (Houde et al. 2002), DR 21 (OH) (Lai, Velusamy, k Langer 
2003), and Orion A (Houde et al. 2004). 

Once a and (3 are mapped across a given molecular cloud, the angular dispersions a (a) 
and a {(5) can be calculated from the measured data. It is easy to show that, in the small 
angle limit, the total angular dispersion of the magnetic field lines a (0) is given by 

a 2 ((f)) = a 2 (a) + sin 2 (a) a 2 ((3) . (5) 

Equation (5) takes into account not only the inclination of the magnetic field, but also 
angular deviations along two independent directions perpendicular to the field orientation. 
Because of this last point, the velocity dispersion will be a/2 times larger than what is used 
in the original CF method equation (4). That is to say, we will now either use the two- 
dimensional velocity dispersion &(v±), define after equation (1), or its equivalent expressed 
as a function of a (v) if the turbulence is isotropic 




Using equations (5) and (6), we can now write a generalized CF equation from (1) 



i 

2 

(7) 



or if the turbulence is isotropic 



B = C 



Aiipa (v±) 
a 2 (a) + sin 2 (a) a 2 W) 
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B = C 



8npa 2 (v) 
3 [a 2 (a) + sin 2 {a) a 2 (/?)] 



(8) 



In both equations (7) and (8) we have added a correction factor C (first introduced 
by Ostriker, Stone, & Gammie (2001)) to take into account some shortcomings of the CF 
method to be discussed later. It is now easy to see how equation (8) can be reduced to one 
for the plane-of-the-sky component of the magnetic field -B pos (i.e., equation (4)) when only 
polarization measurements are available. In this case, for a sufficiently large set of data we 
expect (as long as a ^ 0) 

a 2 (a) = sin 2 (a) a 2 {(5) , 



and 



a 2 (0) = 2 sin 2 (a) a 2 {(3) 



We can write 



Epos = B sin (a) 

fl <t(v) 

which is the same as equation (4). 

We can, therefore, emphasize two important advantages of the modified CF equation 
(7) (or (8)) over the original: 

• the new equation is valid no matter what the orientation of the magnetic field is. Most 
notably, the method does not fail when the field is directed along the line of sight. 

• Finally, the value for the magnetic field calculated with equation (7) is not that of 
its plane-of-the-sky component, but that of full magnitude of the mean magnetic field 
vector. 
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2.2. Shortcomings of the method 

As mentioned earlier, MHD simulations have already been used in the past (Ostriker, 
Stone, k Gammie 2001; Padoan et al. 2001; Heitsch et al. 2001; Kudoh k Basu 2003) to test 
the validity of the original CF method (equation (4)). The main conclusion of these studies 
was that the introduction of a correction factor (C in equations (7) and (8)) is needed when 
evaluating -B pos - A correction of C ~ 0.5 was deemed appropriate in most cases when the 
field is not too weak. A few reasons are usually identified for this. For example: 

1. Smoothing of the field: because of the finite resolution with which observations are 
done, there will be an averaging of the angular structure of the field. This will bring 
a decrease of the angular dispersion a ((f)), and an overestimation of the field strength 
(Ostriker, Stone, k Gammie 2001). 

2. Similarly, line-of-sight averaging (independent of the angular resolution of the obser- 
vations) of the magnetic field will decrease a (0) (Myers k Goodman 1991). 

3. Inhomogeneity and complex density structures (e.g., dumpiness) also tend to reduce 
the value of C (Zweibel 1990; Ostriker, Stone, k Gammie 2001). 

We also add to the previous points one more aspect that should be kept in mind when 
applying the CF method. 

In the case of highly turbulent and massive molecular clouds (like in the example consid- 
ered in the next section), it has been observed that there can exist significant velocity drifts 
between coexistent neutral and ionic molecular species. This can be ascertained through the 
comparison of the observed line profiles for the two types of species, the ions consistently ex- 
hibiting narrower spectral line widths (Houde et al. 2000a, b, 2002; Lai, Velusamy, k Langer 
2003; Houde et al. 2004). This implies that the coupling between ions and neutrals is not 
perfect (Houde et al. 2002). Within the context of the CF method, this bring about uncer- 
tainties in two of the quantities used when evaluating the magnetic field strength. Indeed, 
because of this imperfect coupling between ions and neutrals, the mass density p used in the 
CF equation cannot be that of (larger) neutral density. It must be somewhat smaller. Fur- 
thermore, because of the aforementioned velocity drift, the velocity dispersion perpendicular 
to the field lines cr(v±) (or a (v)) cannot be that measured for a neutral molecular species. 
It must also be smaller. The combination of these factors will also tend to reduce the value 
of C (in equations (7) or (8)), at least when the CF method is applied to highly turbulent 
and massive molecular clouds. 
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We leave the quantification of these effects as open questions that could, perhaps, be 
investigated through simulations. 

3. Application of the extended CF method to M17 

Using their aforementioned technique, Houde et al. (2002) measured the spatial orien- 
tation of the magnetic field at 57 different positions across the M17 molecular clouds. This 
was accomplished using extensive 350 /im dust continuum polarimetry and spectroscopy 
(HCO + /HCN) maps obtained at the Caltech Submillimeter Observatory. We now use their 
results to calculate mean magnetic field strength for M17, using equation (7) 1 . 

From the analysis of Houde et al. (2002) we find the following averages for M17: 

a ~ 47.1°, 
a (a) ~ 10.8°, 

(3 ~ 76.1°, 
a ((3) ~ 16.7° 
a ((f)) ~ 16.3°, 
a(v±) ~ 2.0km/s. 

The transverse velocity dispersion was evaluated from the HCN spectra, taking into 
account the (anisotropic) turbulent flow model used by Houde et al. (2002) (see their Figure 
2, and equation (11)), and the fact that the inclination angle is known 2 . Upon using equation 
(7) with C = 0.5, an approximate value of 10 6 cm -3 for the mean density, and a mean 
molecular mass of 2.3, we find 

B ~ 2.5 mG. 



lr The values for a and (3 used here are slightly different from those presented in Houde et al. (2002). We 
use a maximum polarization level of 10%, instead of 7% as was used in their original analysis. See Houde et 
al. (2004) for more details. 

2 Within the context of the anisotropic turbulent model of Houde et al. (2002), a value for a (v±) at each 
position can be obtained from the corresponding observed spectral line width <r bs{v)- It can be shown 
that a 2 (v±) = a 2 hs (v) f / [ecos 2 (a) + f/2 sin 2 (a)] , where e and / are given in their equation (11) with 
A6» = 44.4°. 
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This value for the magnitude of the magnetic field could be further reduced if the 
correction factor C were found to be smaller than the stated value (because of the effects 
discussed in section 2.2), or again if the average density across the maps were less than 
what was assumed here. However, this field strength may not be too excessive in light of 
the fact that Brogan & Troland (2001) obtained a peak value of — 750 /zG for the line-of- 
sight component of the magnetic field in M17, using HI Zeeman measurements. For, once 
the inclination angle quoted above is taken into account, we calculate from their data a field 
magnitude in excess of 1 mG. Our molecular species (i.e., HCN and HCO + , in the J = 4 — > 3 
transition) probe denser media which could harbor stronger fields. 

It is also interesting to note that 

a (a) ~ sin (a) a {(3) = 12.2°, 

as would be expected for a large enough data set. 

Finally, we would like to state that the extension to the CF method presented in this 
paper should be readily testable through MHD simulations, as was done in the past for the 
original CF technique (Ostriker, Stone, & Gammie 2001; Padoan et al. 2001; Heitsch et al. 
2001; Kudoh & Basu 2003). 

The author thanks T.G. Phillips, R. Peng, and S. Basu for helpful discussions. The 
Caltech Submillimeter Observatory is funded by the NSF through contract AST 99-80846. 
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Fig. 1. — The spatial orientation of the magnetic field is defined with the two angles a and 
(5. The N, E, and LOS axes stand for north, east, and line of sight, respectively. From Houde 
et al. (2002). 



